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In the framework of a (2+l)-dimensional P-even massive Gross-Neveu model, an external magnetic 
field is shown to induce a parity breaking first order phase transition. Possibility of applying the 
results obtained to description of magnetic phase transitions in high-temperature superconductors 
is discussed. 



- . I. INTRODUCTION 

o 

The dynamical symmetry breaking phenomenon induced by external magnetic or chromomagnetic 
fields is called magnetic catalysis effect. This property of a uniform external magnetic field H was 
for the first time observed in the study of (2+l)-dimensional (3-d) chiral invariant field theories 
with four-fermion interaction (the so called Gross - Neveu (GN) type theories Q). In this case, 
an arbitrary weak external magnetic field leads to dynamical chiral symmetry breaking (D^SB) 
even for an arbitrary small coupling constant p],|3[]. This phenomenon was then explained basing 
upon mechanism of effective reduction of space-time dimensions in the external magnetic field and, 
corresponding strengthening of the role of infrared divergences in the vacuum reorganization |3J. 
Later, it has been shown that spontaneous chiral symmetry breaking can be induced by an external 
chromomagnetic field as well |5|-|7| . Moreover, basing upon the study of a number of field theories, it 
has been argued that D^SB magnetic catalysis effect may have a universal, i.e. model independent, 
character || (for 3-d quantum field theories this fact was proved in J3J). In a recent paper |J in the 
framework of a P— even 3-d GN model, it has been demonstrated that the external magnetic field 
serves as a catalyst for spontaneous parity breaking as well. Magnetic catalysis has already found its 
applications in cosmology and astrophysical investigations [|10|| , and also in constructing the theory of 



(N 
c3 



[ high-temperature superconductivity [P],[TTfTr3[| . It can be positively stated that this phenomenon will 
find its application in the elementary-particle physics, condensed-matter physics, physics of neutron 
stars etc., i.e., in those branches of science, where the dynamical symmetry breaking plays a crucial 
role, and an external magnetic field is present. [] For the last ten years, 3-d field theories including GN 
type models |1| have attracted a considerable interest. In many respects, this can be explained by 



the fact that high-temperature superconductivity (HTSC) is a planar phenomenon, i.e., conduction 



electrons in HTSC materials are concentrated in planes formed by atoms of Cu and O ||16|| . In recent 



experimental studies of high temperature superconductor Bi 2 Sr 2 CaCu 2 08 [0, it was discovered that 
thermal conductivity as a function of an external magnetic field H experiences a jump at a certain 
value of H = H c ~ T 2 (the sample temperature T < T c , where T c is the temperature of transition to 



the superconducting state). The authors of |17j assumed that at H = H c a phase transition induced 
by a magnetic field takes place. Phenomenological description of this phase transition (based upon 
the free energy functional of the system) as a first order parity breaking transition was soon proposed 
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flSfl . Moreover, a theoretical (microscopic) explanation of this phenomenon in the framework of two 
3-d models of the GN type based upon the magnetic catalysis effect of the dynamical symmetry 
breaking were presented ||[lT] []. In one of them, a magnetic field induces dynamical breaking of 
chiral symmetry [11], in the other, a parity breaking phenomenon takes place ||. A general feature 



of phase transitions for these two models is their continuity. In other words, both chiral symmetry in 
j 1 | and parity in 0] are violated at the point H = H c by means of the second order phase transition, 



and this does not agree with the phenomenological approach of JT8[ . 

In the present paper, we show that there exists a (2+l)-d GN model, where the magnetic catal- 
ysis demonstrates the same qualitative properties, as in phenomenological description of a phase 
transition in HTSC systems induced by an external magnetic field [[TJ|. 

II. STATEMENT OF THE PROBLEM 

Consider the influence, an external magnetic field exerts on the phase structure of the (2+l)-d GN 
model, described by the Lagrangian 

N - G N 

L = Ys^o.idll) a + rmfj a ^ a ) + ^a^a) 2 , (1) 

a=l a=l 

where d = 7^(9^, the fields ip a are transformed according to the fundamental representation of the 
U(N) group introduced in order to employ the nonperturbative 1/iV-expansion method. Moreover, 
for every value of a = 1,2, ...,N, Fermi fields ip a are four- component Dirac spinors, reducible with 
respect to the Lorentz group (corresponding indices are omitted). The matrix r in the spinor space 
has the form: r = diag(l, 1, — 1, — 1). The action with the Lagrangian (|I]) is invariant under the 
discrete parity (P) transformation i/j a (t, x, y) — > ij j 5 i/) a (t, — x, y) (the algebra of 7^ matrices for the 
reducible 4-dimensional spinor representation of the Lorentz group for 3-d space-time is presented 
in 129.) It is this model at m = that was used in || for the description of phase transitions 
induced by an external magnetic field in HTSC systems. We demonstrate that, in contradistinction 
to the case of m = 0, in this model at m 7^ at the critical point H = H c , the external magnetic 
field induces dynamical breaking of parity that takes place discretely, i.e., in the framework of the 
first order phase transition. This property of the magnetic catalysis in the model |) at m / 
has no analogy in any other field model, and this is demonstarated for the first time in the present 
paper. With the results obtained, one may hope that the massive GN model (P can be used for 
theoretical explanation of the magnetic phase transition in the experiment ||17|| , in agreement with 
the phenomenological approach 



III. PHASE STRUCTURE OF THE MODEL AT H = 

Before considering the influence of the nonzero external magnetic field H on the vacuum of the 
model (HI), its phase structure will be studied at H — 0. To this end, we introduce the auxiliary 
Lagrangian 

- - . - - Na 2 

L = ipidip + mipip + aiprip + , (2) 

2Gr 



2 The relevance of other 3-d relativistic models to this phase transition was recently discussed in 
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where cr(x) is an auxiliary real boson field, and summation over indices a of the auxiliary group 
U(N) is implicitly performed here and in what follows. The field theories ([]]) and (Q) are equivalent, 
since the field a can be excluded from (Q) by means of the equations of motion and the Lagrangian 
(H) is obtained. It is easily shown that the auxiliary field undergoes discrete parity transformation 
P in the following way: a(t,x,y) — ► — a(t, —x,y), i.e. a is a pseudoscalar field. Starting from the 
Lagrangian (0), one can find the effective potential of the theory, which has the following form in 
the one-loop approximation (i.e. in the leading order of the 1/N expansion): 



V (a) = 



Na 2 
2G 



where integration is performed over the Euclidean momentum, and 
(0) over the domain < p 2 < A 2 yields: 



—V (a) = — + 
N 2g 



\m — a\ 
6tt 



+ 



\m + a\ 

67T 



(3) 



m ± a | . Integration in 



(4) 



where: - = ^ — fr- It is well known that the phase structure of any theory is to a great extent 
determined by the symmetry of the global minimum of its effective potential. It is easily shown 
that the function (B) has a global minimum point a = for ir/g > —2m, symmetric with respect to 
parity transformation. If n/g < —2m, potential @) obtains in the domain < a < oo a nontrivial 
absolute minimum at the point 



7T 71- 

2g V 4 ^ 2 



(5) 



and parity of the model is spontaneously broken. When we use a bare coupling constant G insted 
of g, the following picture is obtained. At G < G c = 7r 2 /(2A — 2mn) the vacuum of the model is 
P— even, at G > G c a phase of the model with spontaneously broken parity is realized. At the point 
G = G c the global minimum of the potential jumps from the origin to the nontrivial point (||), and 
the first order phase transition takes place. (We note for comparison that in the case of a massless 
theory, i.e., for m = 0, a second order phase transition continuous in the coupling constant takes 
place at the point G = G c .) 

IV. MAGNETIC CATALYSIS OF DYNAMICAL PARITY BREAKING 



We now consider the influence of an external magnetic field H (at zero temperature) on the 
symmetric phase of the model (0), i.e., in the case n / g > —2m (the bare coupling constant is 
sufficiently small in that case: G < G c ). The corresponding effective potential, which is a special 
case of the effective potential for a more general GN model at H ^ 



has the form: 



V H (a) 



No 2 



+ nj: 



k=l 



eHM k (2eH) 



3/2 



47T 



47T 



-c 



2eH 



(6) 



where ((s,x) is the generalized Riemann zeta-function pi] . Let us demonstrate that an external 
magnetic field induces spontaneous parity breaking for this model as well. Moreover, the magnetic 
catalysis phenomenon and its properties being qualitatively different for cases m = and m^O. 

The case with m = 0. The corresponding effective potential of the model V^ =0 (a) is obtained 
from the formula (|6|) at m = 0: 



3 



Na 2 



NeH\a\ N(2eHfl 2 ( 1 a 



2tt 



2tt 



C 



2' 2eH 



(7) 



It is clear that V% =0 



cr) is a function symmetric under the transformation a — > —a. Hence, in order 
to find a global minimum point, its investigation only in the set a G [0, oo) will suffice. In this case 
the stationary equation 



dV, 



m=0 1 
H 1 



a 



da 



= 2NaF(a) = 2Na 



1 



eH 

47TO" 



V2eH /l 
4vr ^\2'2eH 







(8) 



is obtained from (0) with the help of the formula d((s, x) / dx= — s((s + 1, x). The potential (0) was 
studied in [||,[|||. Nevertheless, we will dwell upon some details that we will need in what follows. 
First of all, we point out that it follows from (Kl) that 



dV, 



m=0i 
H 1 



a 



d. 



a 



= lim (2NaF(a)) 

^0+ a— >0+ 



NeH 



2tt 



(9) 



Moreover, since V, 



1) At the pont a 

rm=0 1 



m=0i 
H 1 



cr 



0, a local 



The first derivative of the function Vff (cr) does not exist 



i.e., the point cr = is not a solution of the stationary equation 
V^ =0 (— cr), there appear two more consequences of the formula ( 
maximum of the potential is situated. 2 
at the point cr = 0. (At the same time, at H = 0, the effective potential is a differentiable function 
on the whole cr— axis.) Consequence 1), as well as the fact that lim^i^cx, V^ 1=0 (cr) = +oo, indicate 
the presence of a nontrivial global minimum of potential V^ =0 (cr) at the point ctq(H) ^ 0. This 
means that parity of the massless model under consideration in the presence of an external magnetic 
field (arbitrary small) is inevitably spontaneously broken even for arbitrary small coupling constant 
G (in this case g > 0). Hence, a dynamical breaking of parity, induced by an external magnetic field 
takes place (magnetic catalysis phenomenon). 

It is known that F(a) is a monotonically increasing function in the set cr 6 [0, oo), so that 
F(0) = — oo (the consequence of equality <^)) and \xm a ^ 00 F(a) = oo Hence, there exists 

a unique point ao(H) ^ 0, where the function F(a) vanishes, and this is the place, where the global 
minimum of potential ([I]) is situated. It follows from the equation F(a) = 0, implicitly defining the 
function <jq(H), that <Jq(H) w egH/(2ir) at H — ► Thus, when a magnetic field is included in 

the model, a parity breaking second order phase transition continuous in H takes place, due to the 
fact that the order parameter cr (H) is a continuous function of H at the point of phase transition 
(i.e., at H = 0). 

It was proved earlier [|2]|| that at large values of H the function a (H) behaves as follows: cr (H) rs 
k\fe~H, where k is a solution of the equation 1 = \/2£(l/2, k 2 /2) and k ~ 0.45. As a (H) is 
the point of global minimum of potential (0), the following inequality is valid for all values of H: 
(0) > V^ =0 '(<7o(H)) . In particular, using this relation for H — ► oo renders 



T/m=0 
" H 



N{2eHf/ 2 (I 
2tt ^ V 2' 



> 



Nk(eH) 
2V 



3/2 



N{2eHf/ 2 (1 If 
2~k ^ V 2' ~2 



(10) 



The case with m^0. We will here discuss some special features of spontaneous symmetry break- 
ing in the model ([!]) at H, m ^ 0. As in the massless case, at m ^ potential Vff(a) is an even 
function of cr, and hence, suffices to study it on the semiaxis cr e [0,oo). Here, the stationary 
equation for the effective potential has the form: 







1 dV, 



H 



N da 



(cr + m)F(a + m) — (m — a)F(m 
(cr + m)F(a + m) + (cr — m)F(a - 



-a) 
m) 



a < m; 
a > m, 



(11) 
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where function F(x) is presented in (|8|). With regard for it is easily obtained from the above 
equation 



dVi 



H 



da 



NeH dV H 
= 1 — 

a^m + 2ir da 



(12) 



This means that at H ^ potential (|6]) is not differentiable at points a = ±m. By virtue of this, 
the possible global minimum point of the function V#(er) in the set a G [0, oo) is either at the point 
a = m, or at one of the solutions of the stationary equation fllTI). However, using (|12"D, one can make 



the important conclusion: if, at a certain vicinity to the left of the point a = m, the derivative of 
the effective potential is negative, its derivative in a certain vicinity to the right of it is also negative. 
Hence, a = m is unable to be not only the global minimum point, but even a local minimum point. 
Therefore, special attention should be paid to solutions of equation ([□]). An evident solution of this 
equation for all values of H is the pont a = 0. By means of both analytical and numerical methods, 
one can demonstrate that for all H ^ there exists at this point at least a local minimum of the 
function Vjj(a). Moreover, for sufficiently small values of H the global minimum of the potential is 
situated just at the point a = 0, and parity remains unbroken (this is the consequence of the fact 
that with such H the derivative of the potential is positive for all a G [0, oo)). 

We will demonstrate that at large H the equation ([□]) has one more solution that is absent for 
small H. Indeed, in the domain a » m, this equation coincides in form with the stationary equation 
(§) for the massless case. The solution cr (H) of the latter becomes large only in the limit eH — > oo. 
Hence, in the region of large values of the magnetic field, equation (|TTD has, besides a = 0, one more 
solution ao(H), such that 5q(H) ~ kyfeH at eH — > oo (the value of coefficient k is the same as in 
([10|)). Besides a = and ao(H), we were able to find no more solutions of the stationary equation 

As it was already pointed out above, at sufficiently small H, the global minimum of the potential 
is situated at the point a = 0. We will demonstrate that with growing H it goes over to the point 
a (H). To this end, we will find the values of the potential in stationary points for eH/m 2 — > oo: 

Vh(0) = - (2e f )3/2 C (-1 O) + 0(eH/m 2 ), (13) 



2tt a V 2 

Vh(MH)} = « _ (M^ C (_!,£) +CW , (14) 



where we considered that &q(H) ~ kv eH for eH — > oo. Comparing expressions ( |T3| ) and ( |HD with 
the help of fllUp, we come to the conclusion that Vh{vq(H)) < Vh(0) at eH — > oo. Since for all 
values of H, a = is at least a local minimum, the passing of a global minimum of the potential 
from the point a = to the nontrivial point (Tq{H) takes place in a jump at a certain critical value 
of an external magnetic field H c (g) ^ 0. Thus, at the point H c (g) a parity breaking first order phase 
transition takes place. The critical magnetic field H c (g) as a function of the coupling g at T = and 
m 7^ is depicted in Fig. 1. In Fig. 2, the functions Vh(c) are plotted at various values of H and 
gm = 10. 

V. MAGNETIC CATALYSIS AT NONZERO TEMPERATURE 

Employing the technique developed in || , one can obtain the following expression for the effective 
potential of the model fll]) at H, T ^ 0: 



5 



1 1 g TT 2 OO 

jjV HT (a) = -V H {a) - — £{m(l + exp(-M^)) + 2 £ ln(l + exp(-(3^Mf + 2keH))}, (15) 

where j3 = 1/T, and Vh(g) is the potential Numerical analysis of the potential ( JT5| ) demonstrates 
that also in this case parity is unbroken at sufficiently small values of H. However, there exists a 
critical magnetic field, when a parity breaking first order phase transition appears. At fixed value 
of g, we denote the critical field value as H C (T). Some of its values for gm = 5 are presented in the 
Table 1. From this Table it is seen that H C (T) ~ T 2 at T — > oo, i.e., H C (T) behaves as in models 
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VI. CONCLUSIONS 

In the present paper, it was proved, in the framework of quantum field theory, that an external 
magnetic field can induce a first order phase transition. As an illustration, we proposed a P— even 
massive 3-d GN model (|l|) |. 

At m = 0, magnetic catalysis of the dynamical parity breaking takes place in this model 0. 
Moreover, at the point H c ~ T 2 (T — ► oo), there appears a second order phase transition, and 
theoretical values i? c provide satisfactory description of the experimental data [|T7jl . 

We have shown that, in the model ([1]) at m ^ 0, dynamical parity breaking is also induced by an 
external magnetic field. The curve H c at large T behaves in the same way as at m = 0. However, in 
contradistinction to the latter case, in the massive GN model ([!]) at the point H c , there appears a 
parity breaking first order phase transition, and this fact is in agreement with the phenomenological 
description ]18[ of magnetic phase transitions in HTSC systems in the experiment |17| . 
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FIG. 1. Critical magnetic field as a function of coupling constant g 
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FIG. 2. Effective potential at various values of magnetic field intensity and gm = 10 (a/m — x). 



TABLE I. Some values of H C {T) at gm = 5. 



T/m 


0.1 


0.25 


0.5 


1 


2.5 


5 


10 


25 


50 


100 


200 


eH c (T)/m z 


24.9 


25.2 


27 


36 


131 


465 


1782 


10900 


43200 


173000 


700000 
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